Stochastic Processes Crossing from Ballistic to Fractional Diffusion with Memory: 

Exact Results 
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We address the now classical problem of a diffusion process that crosses over from a ballistic 
behavior at short times to a fractional diffusion (sub- or super-diffusion) at longer times. Using the 
standard non-Markovian diffusion equation we demonstrate how to choose the memory kernel to 
exactly respect the two different asymptotics of the diffusion process. Having done so we solve for the 
probability distribution function (pdf) as a continuous function which evolves inside a ballistically 
expanding domain. This general solution agrees for long times with the pdf obtained within the 
continuous random walk approach but it is much superior to this solution at shorter times where 
the effect of the ballistic regime is crucial. 



Introduction: Nature offers us a large number of ex- 
amples of diffusion processes for which an observable X 
diffuses in time such that its variance grows according to 
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for t <C i c j 
for t^>t c . 
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where angular brackets mean an average over repeated 
experiments and D 2 and D a are coefficients with the 
appropriate dimensionality. The short time behavior is 
known as 'ballistic', and is generic for a wide class of 
processes. The long time behavior with a ^ 1 is generic 
when the diffusion steps are correlated, with persistence 
for a > 1 and anti-persistence for a < 1 [l| . These corre- 
lations mean that the diffusion process is not Markovian, 
but rather has memory. Thus the probability distribution 
function (pdf) of the observable X, f(X,t) is expected 
to satisfy a diffusion equation with memory 0, 



d.f(X,t) 
dt 



dt'K(t-t')X7 2 f(X,t) 



(3) 



with Kit) being the memory kernel and V 2 the Laplace 
operator. 

In this Letter we study the class of processes which 
satisfy Eqs. ([TJ>- (J3J) - First of all we find an expression 
for the kernel K(t) which is unique for a given law of 
mean-square-displacement. Second we consider the ker- 
nel which contains both the ballistic contribution embod- 
ied in Eq. ([1]) and the long-time behavior ((2]). For this 
case we find an exact equation and a solution for Eq. 
(|3|). Lastly a simple interpolation formula for the kernel 
is inserted to the exact equation which is then solved for 
the pdf of X without any need for the fractional dynam- 
ics approach [jj. Some interesting characteristics of the 
solution are described below. 

Determination of the kernel K{t): To determine 
the kernel in Eq. §S§ we use a result obtained in 
Consider the auxiliary equation 



dP(X,t) 
dt 



V 2 P{X,t) 



(4) 



Define the Laplace transform of the solution of Eq. ((4]) 

as 



P(X,s) 



dte~ st P(X,t) , 



(5) 



it was shown in Q that the solution of Eq. with the 
same initial conditions can be written as 



f(X,s) 



1 p (x 3 
K(s) ' K(s) 



), 
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where here and below the tilde above the symbol means 
the Laplace transform. The development that we pro- 
pose here is to replace in Eq. ([6|) the Laplace transform 
K(s) with the Laplace transform of the mean-square dis- 
placement. This is done by first realizing (by computing 
the variance and integrating by parts) that 



d(X 2 )(t) 



dt 



2 [ dt'K(t-t') 



or, equivalently, 



K(s) 
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The second line was written in order to find the time 
representation of K(t) which is the inverse Laplace trans- 
form: 

*(*) = K^I <X2>(<) + £ <X2>W ) 



ld_ 

2dt 



H(t) 



d(X 2 )(t) 
dt 



(9) 



where H(t) is the Heaviside function. Obviously, using 
the first line of Eq. ([8]) in Eq. ([6]) the solution is entirely 
determined by whatever law is given for the variance, 
together with initial conditions. 
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For ordinary diffusion the variance is defined by Eq. 
© with a = 1 and t c = 0. It follows from Eq. © 
that the kernel is K(t) ~ and Eq. ([3]) is reduced 
to the Markovian Eq. ©; this process does not possess 
any memory. More complicated examples are considered 
below. 

Example I: fractional differential equations. In 

recent literature the problem of a diffusion process which 
is consistent with Eq. © only for all times (i.e. t c = 0) is 
investigated using the formalism of fractional differential 
equations (see, e.g., (3]). In this formalism Eq. © is 
replaced by the fractional equation 



df(X,t) 
dt 



- a d 2 f(X,t) 
dx 2 



(10) 



where the Rieman-Liouville operator qDj a is defined by 



T(a)dtJ (t-t') 1 -"' 



(11) 



where T(a) is the gamma function. It is easy to see 
that this equation follows from Eq. © with the kernel 
evaluated by Eq. © with the variance ©. We reiterate 
however that this equation is consistent with Eq. © 
for all times t > 0. This of course is a problem since 
this formalism cannot agree with the ballistic short time 
behavior which is generic in many systems. 

Example II: ballistic behavior. For X one- 
dimensional the solution of Eq. © with the initial con- 
dition P(X, t = 0) = S(X) is given by 



P(X,s) 



2^ 



exp(- | X | Vi) 



Substituting in Eq. © we find 



f(X,s) 



1 



2s 3 (X 2 )(s) 



: exp(- | X 



s(X*)( S y 



(12) 
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For systems with the pure ballistic behavior (e.g., dilute 
gas) the variance can be written as (X 2 ) t = (u 2 )t 2 , where 
(u 2 ) is the mean square average of the particle veloci- 
ties. The Laplace transform of this expression is given 
by (X 2 )(s) = 2(u 2 ) / s 3 and the Laplace transform of the 
pdf is defined by 



f(X,s) 



1 



exp( 



X 



vV>' 



(14) 



The inverse transform reads 



f(X,t) = ls(\X\-^)t) 



(15) 



This solution corresponds to a deterministic evolution; 
there is a complete memory of the initial conditions in 
the absence of inter-particle interactions, Kit) = (u 2 )). 

General case: In the general case the mean-square- 
displacement satisfied some law (X 2 )(t) which is sup- 
posed to be known at all times, with possible asymptotic 
behavior as shown in Eqs. © and©. To find the ap- 
propriate general solution we will split f(X, s) into two 
parts, fi(X,s) and fu(X,s), such that the first part 
is constructed to agree with the existence of a ballistic 
regime. Suppose that in that regime, at short time, the 
mean-square-displacement can be expanded in a Taylor 
series 



(X 2 )(t) =J2^ t ^ 1 = a Q t 2 + ai t 3 + a 2 t 4 



(16) 



where /io = 3, /ii = 4 etc. Then the Laplace transform 
(X 2 )(s) can be written for s — > oo as Q 
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Substituting Eq. dTTJ) up to 0(s" 4 ) in Eq. (13]) yields 



fi(X,s) 



1 



exp{- 



X 



is 



3ai 



2,/ao ^/ao 2ao 

The inverse Laplace transform of this result reads 



)) • (18) 



fj(X,t) = \ex P ^t)8i\ X 
i Ado 



dot) 



(19) 



Not surprisingly, this partial solution corresponds to a 
deterministic propagation. Note that in order to avoid 
exponential divergence in time we must have a\ < in 
the expansion (|T6|) . 

Having found fiiX, s) we can now write ///(X, s) sim- 
ply as 



f II (X,s)=fiX,s)-f I iX,s) . 
Calculating this difference explicitly we find 



(20) 



fniX,s) = -I 



s 3 (X 2 )is)) 



exp(- \X\{ 



<X 2 )( 6 



))- 



1 . 3ai 
exp( 



2a 



3/2 



X 



cxp( _L^ s ) 



FiX, s) exp(- 



X 



s). 



(21) 
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The inverse Laplace transform of Eq. (|21[) is given by 
f I i(X,t)=F(x,t-\X\/y/ES)H( y /ESt-\X\). (22) 



The importance of this result is that the explicit Heav- 
iside function is taking upon itself the discontinuity in 
the solution fjj(X,t). The exact value of this function 
at the point | X \= ^/ao~t can be calculated using the 
initial value theorem and is given by 



///(I X | 
/ 3 ai 
4' 



a t, t) 



1 ,27 r a x . 



3/2 
V 



2Joq~ 8 a 




(23) 

x exp(--— t) . 
I a 



Summing together the results (JT9J) and (|22|) in the time 
domain we get a general solution of the non-Markovian 
problem with a short-time ballistic behavior, in the form 



f(X,t) 



F(X,t- 



X 



)H(yfot- | X |). (24) 



This is the main result of the present Letter. The diffu- 
sion repartition of the probability distribution function 
occurs inside the spatial diffusion domain which increases 
in a deterministic way. The first term in Eq. (j2"4"|) corre- 
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FIG. 1: The time evolution of the function fi(X, t) defined by 
Eq. (|19p for time intervals t/to —0.5, 1, 2, 4, 8 (the time scale 
to — ao/(3ai)). The S- function is graphically represented by 
narrow Gaussians. 

sponds to the propagating <5-function which is inherited 
from the initial conditions, and it lives at the edge of the 
ballistically expanding domain. Schematically the time 
evolution of this term is shown in Fig[TJ where the 6- 
function is graphically represented as a narrow Gaussian. 
The dashed line represents the exponential decay of the 
integral over the (5-function. The function F(X, t — ^=) 
in the time domain is a continuous function and can be 
evaluated numerically, for example using the direct inte- 
gration method 0. Below we demonstrate this calcula- 
tion with explicit examples. 



Interpolation for all times: To interpolate Eqs. (Q]) 
and ([2]) we propose the form 



(25) 



where < a < 2. Here to is the crossover characteristic 
time, at t <§; to the law P5|) describes the ballistic regime 
and at t ^> to the fractional diffusion. 





FIG. 2: The continuous part of the pdf (|22[) for different values 
of the parameter a. Superdiffusion (a = 3/2, upper panel), 
regular diffusion (a = 1, middle panel) and subdiffusion (a — 
1/2, lower panel). Time intervals from the top to the bottom 
r =0.5, 1, 2, 4, 8. The reader should note that the full solution 
of the problem is the sum of the two solutions shown in this 
and the previous figure. 

Introduce now dimensionless variables (£ 2 )t = 
(AX 2 ) t /(2D a t^) and r = t/t . With these variables 
the last equation reads 



(l + r) 



2-Q ' 



(26) 
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The Taylor expansion of (j25|) is given by 
(£ 2 >r = r 2 - (2 - a)r 3 + 1(3 - a) (2 - a)r 4 + . . . . (27) 

Substitution these expansion coefficients into Eq. (fl~9"l) 
yields the first term in the expression for the probability 
distribution function ([24| 

f r (x,t) = lexp{-^^-r)5{\ £ \ -r). (28) 
The Laplace transform of Eq. (|2"rJ)) is 

s s \ s J s a 

(29) 

where T(a, s) is the incomplete gamma function. Note 
that the case a = 2 is special, since it annuls the ex- 
ponent in Eq. (|28|) . leaving as a solution a ballistically 
propagating 5-function. For all other values of a < 2 the 
inverse Laplace transform of the function F(x, s) which 
defines the diffusion process inside the expanded spatial 
domain should be evaluated, in general, numerically. 

Results of the calculations following the method of Ref. 
6] for the smooth part of the probability distribution 
function fu(x,t) for different values of the parameter 
a are shown in Fig. [5J The reader should appreciate 
the tremendous role of memory, or the non-Markovian 
nature of the process under study. For example regu- 
lar diffusion with a = 1 results in a Gaussian pdf that 
is peacefully expanding and flattening as time increases. 



Here, in the mid panel of Fig. [2] we see that the ballistic 
part which is represented by the advancing and reduc- 
ing (5-function sends backwards the probability that it 
loses due to the exponential decay seen in Fig. [T] This 
'back-diffusion' leads initially to a qualitatively different 
looking pdf, with a maximum at the edge of the ballis- 
tically expanding domain. At later times the pdf begins 
to resemble more regular diffusion. The effect strongly 
depends on a simply due to the appearance of a in the 
exponent in Eq. (|28|) . 

For long times the solutions shown in Fig. [2] agree with 
the Markovian pdf obtained in the frame of a continuous- 
time random walk Q. For the special case a = the 
limiting behavior of the general solution from Eq. (j2"4"|) 
can be evaluated with the help of the final value theorem: 

f(X) = ^=eM-V2\X\), (30) 

this analytical result coincides with the pdf from [7fl at 
the same conditions. 

In summary, we have shown how to deal with diffusion 
processes that cross-over from a ballistic to a fractional 
behavior for short and long times respectively, within 
the time non-local approach. The general solution (f24|) 
demonstrates the effect of the temporal memory in the 
form of a partition of the probability distribution func- 
tion inside a spatial domain which increases in a deter- 
ministic way. The approach provides a solution that is 
valid at all times, and in particular is free from the in- 
stantaneous action puzzle. 
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